SISSA 78/2005/EP 
IC/2005/102 
TUM-HEP-606/ 05 
|hep-ph/0510404| 



The See-Saw Mechanism, Neutrino Yukawa Couplings, 
LFV Decays U — > lj + 7 and Leptogenesis 



S. T. Petcov a ' b >, W. Rodejohann c ), T. Shindou"^ 1 ", Y. Takanishi d )* 

a ) Scuola Internazionale Superiore di Studi Avanzati, 1-34014 Trieste, Italy 

fe ) Istituto Nazionale di Fisica Nucleare, Sezione di Trieste, 1-34014 Trieste, Italy 

c ) Physik-Department, Technische Universitdt Miinchen, 
James-Franck-Strafie, D-85748 Garching, Germany 

rf ) The Abdus Salam International Centre for Theoretical Physics, 
Strada Costiera 11, 1-34100 Trieste, Italy 



Abstract 

The LFV charged lepton decays /1 — > e+7, t — ► e+7 and t — > /i+7 and thermal leptogenesis 
are analysed in the MSSM with see-saw mechanism of neutrino mass generation and soft SUSY 
breaking with universal boundary conditions. The case of hierarchical heavy Majorana neutrino 
mass spectrum, M\ <C M 2 <C M 3 , is investigated. Leptogenesis requires M\ ^ 10 9 GeV. 
Considering the natural range of values of the heaviest right-handed Majorana neutrino mass, 
M 3 >5x 10 13 GeV, and assuming that the soft SUSY breaking universal gaugino and/or scalar 
masses have values in the range of few x 100 GeV, we derive the combined constraints, which 
the existing stringent upper limit on the — > e + 7 decay rate and the requirement of successful 
thermal leptogenesis impose on the neutrino Yukawa couplings, heavy Majorana neutrino masses 
and SUSY parameters. Results for the three possible types of light neutrino mass spectrum - 
normal and inverted hierarchical and quasi-degenerate - are obtained. 
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1 Introduction 



The experiments with solar, atmospheric, reactor and accelerator neutrinos ^ HI El EJ 0>] have 
provided during the last several years compelling evidence for the existence of non-trivial 3-neutrino 
mixing in the weak charged-lepton current (see, e.g., 0): 

3 

viL = ^2UijV jL , l = e,fi,r, (1) 

3=1 

where fn are the flavour neutrino fields, fji, is the field of neutrino Uj having a mass nij and U 
is the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix [Jj, U = J/pmns- The existing 
data, including the data from the 3 H /3-decay experiments 8 imply that the massive neutrinos 
Vj are significantly lighter than the charged leptons and quarks: raj < 2.3 eV (95% C.L.) . A 
natural explanation of the smallness of neutrino masses is provided by the see-saw mechanism of 
neutrino mass generation JTI] . The see-saw mechanism predicts the light massive neutrinos Uj to be 
Majorana particles. An integral part of the mechanism are the heavy right-handed (RH) Majorana 
neutrinos |12| . In grand unified theories (GUT) the masses of the heavy RH Majorana neutrinos 
are typically by a few to several orders of magnitude smaller than the scale of unification of the 
electroweak and strong interactions, Mqut — 2 x 10 16 GeV. In this case the CP-violating decays of 
the heavy RH Majorana neutrinos in the Early Universe could generate, through the leptogenesis 
scenario, the observed baryon asymmetry of the Universe |13| . 

The existence of the flavour neutrino mixing, eq. lfl|). implies that the individual lepton charges, 
Li, I = e,/j,,T, are not conserved (see, e.g., ^U), and processes like pT — ► e~+7, \x~ — > e~ + e + + e~, 
r~ — > e~ + 7, r~ — > \i~ + 7, pT + (A, Z) — ► e~ + (A,Z), etc. should take place. Stringent 
experimental upper limits on the branching ratios and relative cross sections of the indicated 
|AL|| = 1 decays and reactions have been obtained ^3E3 El (90% C.L.): 

B(/i -> e + 7) < 1.2 x HT 11 , B(/i -> 3e) < 1.2 x 1(T 12 , R(/i~ + Ti -> e~ + Ti) < 4.3 x 10" 12 , 
B(t -» At + 7) < 6.8 x 1(T 8 , B(r -» e + 7) < 1.1 x 10" 7 . 

(2) 

Future experiments with increased sensitivity can reduce the current bounds on B(/x — ► e + 7), 
B(r — > /i + 7) and on R(/i~ + (A, Z) — > e _ + (A, Z)) by a few orders of magnitude (see, e.g., [T^]). 
In the experiment MEG under preparation at PSI |19j it is planned to reach a sensitivity to 

B(/i -» e + 7) ~ (10^ 13 - 10~ 14 ) . (3) 

It has been noticed a long time ago that in SUSY (GUT) theories with see-saw mechanism 
of neutrino mass generation, the rates and cross sections of the LFV processes can be strongly 
enhanced |20j . If the SUSY breaking occurs via soft terms with universal boundary conditions at 
a scale Mx above the RH Majorana neutrino mass scale Mr, Mx > Mr 2 , the renormalisation 
group (RG) effects transmit the LFV from the neutrino mixing at Mx to the effective mass terms 
of the scalar leptons at Mr even if the soft SUSY breaking terms at Mx are flavour symmetric 
and conserve the lepton charges L\. As a consequence of these RG-induced new LFV terms in the 

1 More stringent upper limit on rrij follows from the constraints on the sum of neutrino masses obtained from 
cosmological/astrophysical observations, namely, the CMB data of the WMAP experiment combined with data from 
large scale structure surveys (2dFGRS, SDSS) 0: J2j m i < (°- 7 ~ 2 -°) eV ( 95 % C.L.), where we have included a 
conservative estimate of the uncertainty in the upper limit (see, e.g., \W]). 

2 The possibility of "flavour-blind" SUSY breaking of interest is realised, e.g., in gravity-mediation SUSY breaking 
scenarios (see, e.g., [21]). 
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effective Lagrangian at Mr < Mj, the LFV processes can proceed with rates and cross sections 
which are within the sensitivity of presently operating and future planned experiments |20l I22| 
(see also, e.g., [2S1 EH ESI EH1 EH EH1 EH1 EH| ) • In contrast, in the non-supersymmetric case, the 
rates and cross sections of the LFV processes are suppressed by the factor [SJJ (see also [SI]) 
{m j /M w ) i < 6.7 x 10~ 43 , M w being the W ± mass, which renders them unobservable. 

One of the basic ingredients of the see-saw mechanism is the matrix of neutrino Yukawa cou- 
plings, Y u . Leptogenesis depends on Y u as well pS] (see also [SSI CH| and the references quoted 
therein) . In the large class of SUSY models with see-saw mechanism and SUSY breaking mediated 
by flavour-universal soft terms at a scale Mx > Mr we will consider, the probabilities of LFV 
processes also depend strongly on Y u (see, e.g., [2S1 El])- The matrix Y„ can be expressed in 
terms of the light neutrino and heavy RH neutrino masses, the neutrino mixing matrix l^pmnS; 
and an orthogonal matrix R [2S]- Leptogenesis can take place only if R is complex. The ma- 
trix Y v depends, in particular, on the Majorana CP-violation (CPV) phases in the PMNS matrix 
t^PMNS [SH] 3 - It was shown in [201 ISO] that if the heavy Majorana neutrinos are quasi-degenerate 
in mass, the Majorana phases can affect significantly the predictions for the rates of LFV decays 
fx — > e + 7, t — > e + 7, etc. in the class of SUSY theories of interest. 

The matrix Y v can be defined, strictly speaking, only at scales not smaller than Mr. The 
probabilities of LFV processes depend on Yj, cit the scale -^ij? — (Mr). In order to evaluate 
Y u (Mr) one has to know, in general, the light neutrino masses rrij and the mixing matrix £/pmns 
at Mr, i.e., one has to take into account the renormalisation group (RG) "running" of rrij and 
f/pMNS from the scale Mz ~ 100 GeV, at which the neutrino mixing parameters are measured, to 
the scale Mr (see, e.g., [401 130j and the references quoted therein). However, if the RG running of 
uij and £7pMNS is sufficiently small, Y u (Mr) will depend on the values of the light neutrino masses 
rrij and the mixing angles and CP-violation phases in C/pmns a t the scale Mz- 

Working in the framework of the class of SUSY theories with see-saw mechanism and soft SUSY 
breaking with flavour- universal boundary conditions at a scale Mx > Mr, we investigate in the 
present article the combined constraints, which the existing stringent upper limit of the \x — > e + 7 
decay rate and the requirement of successful thermal leptogenesis impose on the neutrino Yukawa 
couplings, heavy Majorana neutrino masses and on the SUSY parameters. The case of hierarchical 
heavy Majorana neutrino mass spectrum, M\ <C M2 <C Ms, is considered. Leptogenesis requires 
Mi <; 10 9 GeV. The analysis is performed assuming that the heaviest RH Majorana neutrino has a 
mass M3 ^5x 10 13 GeV, and that the soft SUSY breaking universal gaugino and/or scalar masses 
(at the scale Mx) have values in the range of few x 100 GeV. One typically gets M3 ^ 5 x 10 13 
GeV in SUSY GUT theories with see-saw mechanism of neutrino mass generation (see, e.g., |41j). 
If the SUSY breaking universal gaugino and/or scalar masses have values in the few x 100 GeV 
range, supersymmetric particles will be observable in the experiments under preparation at the 
LHC (see, e.g., 02]). We find that under the indicated assumptions, the existing stringent upper 
limit on the [i — > e + 7 decay rate cannot be satisfied, unless the terms proportional to M3 in the 
11 — ► e + 7 decay amplitude are absent or strongly suppressed. The requisite "decoupling" of the 
terms oc M3 from the 11 — > e + 7 decay amplitude is realised if the matrix R has a specific form 
which admits a parametrisation with just one complex angle. Using the latter we obtain results 
for the three types of light neutrino mass spectrum - normal and inverted hierarchical (NH and 

3 Obtaining information about the Majorana CPV phases in the PMNS matrix !7pmns if the massive neutrinos are 
proved to be Majorana particles would be a remarkably challenging problem. The oscillations of flavour neutrinos, 
vi — > vy and vi —+ vy , 1,1' — e,(i,T, are insensitive to the two Majorana phases in (7pmns |35l 1361 . The only feasible 
experiments that at present have the potential of establishing the Majorana nature of light neutrinos Vj and of 
providing information on the Majorana phases in (7pmns are the experiments searching for neutrinoless double beta 
((/3/?V)-decay, (A, Z) (A,Z + 2) + e~ + e~ (see, e.g., [Til I571 1551 133 ] ). 
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IH), and quasi-degenerate (QD). For each of the three types of spectrum we derive the leptogenesis 
lower bound on the mass of the lightest RH Majorana neutrino M\. The lower bounds thus found 
in the cases of IH and QD spectrum are ~ 10 13 GeV. The upper limit on B(/j, — > e + 7) in these two 
cases can be satisfied for specific ranges of values of the soft SUSY breaking parameters implying 
relatively large masses of the supersymmetric particles. Using these soft SUSY breaking parameters 
we derive predictions for B(/i — > e + 7), B(r — > e + 7) and B(r — > /i + 7) which are compatible with 
the requirement of successful leptogenesis. 

Our analysis is performed under the condition of negligible RG effects for the light neutrino 
masses rrij and the mixing angles and CP-violation phases in f/pMNS • The RG effects in question are 
negligible in the class of SUSY theories we are considering in the case of hierarchical light neutrino 
mass spectrum (see, e.g., [1U1 EH]). The same is valid for quasi-degenerate Vj mass spectrum 
provided the parameter tan (5 < 10, tan (3 being the ratio of the vacuum expectation values of the 
up- and down-type Higgs doublet fields in SUSY extensions of the Standard Theory. 

2 General Considerations 

2.1 Neutrino Mixing Parameters from Neutrino Oscillation Data 

We will use the standard parametrisation of the PMNS matrix C/pmns (see, e.g 

( C12C13 S12C13 si3e" i<5 \ 



f^PMNS 



-•S12C23 - Ci2S23Si3e lS C12C23 - Si2S23Sl3e lS S23C13 
\ S12S23 ~ Ci 2 C23Sl3e lS -C12S23 - Sl2C23Sl3e i<5 C23C13 / 



+M 



diag(l,e i 2, e i 2 ), (4) 



where Cjj = cos 9ij, Sij = sin 9ij, the angles 9{j = [0, vr/2], 5 = [0, 2ir] is the Dirac CP-violating 
phase and a and (5m are two Majorana CP-violation phases |35| 143] . One can identify the neutrino 
mass squared difference responsible for solar neutrino oscillations, Attt-q, with Am^i = — ml, 
Atti-q = Am^ > 0. The neutrino mass squared difference driving the dominant — > v T (p^ — > v T ) 
oscillations of atmospheric (p^) is then given by \Am 2 A \ = |Am^| = |Am| 2 | > Am^. The 
corresponding solar and atmospheric neutrino mixing angles, #0 and 9 a , coincide with 9\2 and #23 > 
respectively. The angle #13 is limited by the data from the CHOOZ and Palo Verde experiments |44| . 

The existing neutrino oscillation data allow us to determine Am^ , | Am^ | , sin 2 9\2 and sin 2 2#23 
with a relatively good precision and to obtain rather stringent limits on sin 2 #13 (see, e.g., [2*1 1451 US] ^. 
The best fit values of Am^j, sin 2 9\2, (Am^j and sin 2 2^23 read 4 : 

Am^ = 8.0 x 10- 5 eV 2 , sin 2 9 12 = 0.31 , (5) 



I Amfi I = 2.1 x 10" 3 eV 2 , sin 2 26> 23 = 1.0 , (6) 
A combined 3-v oscillation analysis of the solar neutrino, KamLAND and CHOOZ data gives |45j 

sin 2 ^13 < 0.024 (0.044), at 95% (99.73%) C.L. (7) 

The neutrino oscillation parameters Am^, sin 2 #12, (Am^l and sin 2 2^23 are determined by the 
existing data at 3<r with an error of approximately 12%, 24%, 50% and 16%, respectively. These 
parameters can (and very likely will) be measured with much higher accuracy in the future (see, 

4 The data imply, in particular, that maximal solar neutrino mixing is ruled out at ~ 6a; at 95% C.L. one finds 
cos 200 > 0.26 45 , which has important implications |47|. 
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e.g., 6 ). In all further numerical estimates we use the best fit values of Am^, sin 2 #12, (Am^l 
and sin 2 2^23- Whenever the parameter sin #13 is also relevant in the calculations, we specify the 
value used. 

The sign of Am^ = Am^, as it is well known, cannot be determined from the present (SK 
atmospheric neutrino and K2K) data. The two possibilities, Am 2 1 ^ 32 ^ > or Am^^j < corre- 
spond to two different types of z/-mass spectrum: 

- with normal hierarchy mi < 1112 < rn.3, Am^ = Am^ > 0, and 

- with inverted hierarchy m.3 < mi < 7712, Am| = Am 2 2 < 0. 

Depending on the sign of Am\, sgn(Am 2 , L ), and the value of the lightest neutrino mass, min(?rij), 
the u-mass spectrum can be 

- Normal Hierarchical: mi <C m2 <C m%, m<i = (Am|,)2 ~ 0.009 eV, m.3 = |Am^| 5 ~ 0.045 eV; 

- Inverted Hierarchical: m 3 <C mi < TO2, with mi 2 ^ \Am\\ 2 ~ 0.045 eV; 

- Quasi- Degenerate (QD): mi = m2 — m^ = m, m 2 ^> |Am^J, m <; 0.10 eV. 

2.2 The See-Saw Mechanism and Neutrino Yukawa Couplings 

We consider the minimal supersymmetric standard model with RH neutrinos and see-saw mech- 
anism of neutrino mass generation (MSSMRN). In the framework of MSSMRN one can always 
choose a basis in which both the matrix of charged lepton Yukawa couplings, Ye, and the Ma- 
jorana mass matrix of the heavy RH neutrinos, Mjj, are real and diagonal. Henceforth, we will 
work in that basis and will denote by Dn the corresponding diagonal RH neutrino mass matrix, 
D N = diag(Mi, M 2 ,M 3 ), with Mj > and Mi < M 2 < M 3 . The largest mass M 3 will be stan- 
dardly assumed to be of the order of, or smaller than, the GUT scale Mqut — 2 x 10 16 GeV. 

Below the see-saw scale, Mr = min(Mj), the heavy RH neutrino fields Nj are integrated out, 
and as a result of the electroweak symmetry breaking, the left-handed (LH) flavour neutrinos 
acquire a Major ana mass term: 

C?m = -\v% K) ife VLk + h.c. , (8) 

where u^- = C(pLj) T and 

(m v f = vl {Y?y\M- x f\Y v f . (9) 

Here v u = vsinfi, where v = 174 GeV and tan/3 is the ratio of the vacuum expectation values 
of up-type and down-type Higgs fields, and Y u is the matrix of neutrino Yukawa couplings. The 
neutrino mass matrix m„ is related to the light neutrino masses mj and the PMNS mixing matrix 
as follows 

(m u yi = (U*y k m k (U^) kj . (10) 

Using Q and ()10j) . we can rewrite the "matching condition" at the energy scale Mr in the form 

U* U f = v 2 u YjM^Y, . (11) 

where Dj, = diag(mi, m 2 , ffl.3). Thus, in the basis in which the RH neutrino mass matrix is diagonal 
Mn = Dn, the matrix of neutrino Yukawa couplings at Mr can be parametrised as |23| 

Y V (M R ) = — R . (12) 
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Here R is a complex orthogonal matrix 5 R r R = 1. 

In what follows we will investigate the case when the RG running of rrij and of the parameters 
in t/pMNS from Mz to Mr is relatively small and can be neglected. This possibility is realised in 
the class of theories under discussion for sufficiently small values of tan/? and/or of the lightest 
neutrino mass mm(rrij) (HOI; e.g., for tan/3 ^ 10 and/or min(mj) ^ 0.05 eV. Under the indicated 
condition T) v and U in eq. 1)12(1 can be taken at the scale ~ Mz, at which the neutrino mixing 
parameters are measured. 

As is well-known and we shall discuss further, in the case of soft SUSY breaking mediated by 
soft flavour-universal terms at Mx > Mr, the predicted rates of LFV processes such as \i — ► e + 7 
decay are very sensitive to the off-diagonal elements of 

Yt(M R )Y u (M R ) = \ \J^JW V Rt B N R y^Ut , (13) 

v 

while leptogenesis depends on (see also [SSI EU arid the references quoted therein) 

Y„ (Mr) Yt (Mr) = \ R D, R f V 7 ^ • (14) 

In such a way, the matrix of neutrino Yukawa couplings Y y connects in the see-saw theories the 
light neutrino mass generation with leptogenesis; in SUSY theories with SUSY breaking mediated 
by soft flavour- universal terms in the Lagrangian at Mx > Mr, ~Y u links the light neutrino mass 
generation and leptogenesis with LFV processes (see, e.g., 127] ). 



2.3 The LFV Decays k ->• lj + 7 

As was indicated in the Introduction, in the class of theories we consider, one of the effects of RG 
running from Mx to Mr < Mx is the generation of new contributions in the amplitudes of the LFV 
processes [213 122] ■ In the "mass insertion" and leading-log approximations (see, e.g., [221 12U 128 ), 
the branching ratio of li — ► h + 7 decay due to the new contributions has the following form 



B(k -> lj + 7) 



evv) 



total 



(k) 



G 2 F m% 



(3 + al)ml 



8tt 2 



E 



In 



Mx 
M k 



(Y„) 



tan 2 /3 , (15) 



where i 7^ j = 1,2,3, h,h,h — e ,l Jl -, T > m o an( i -^-o = a o m o are the universal scalar masses and 
trilinear scalar couplings at Mx and ms represents SUSY particle mass. It was shown in [2HI that 
in most of the relevant soft SUSY breaking parameter space, the expression 



m 8 s ~ 0.5 rag m\i2 (m,Q + 0.6 my 2 Y 



(16) 



mi/2 being the universal gaugino mass at Mx, gives an excellent approximation to the results 
obtained in a full renormalisation group analysis, i.e., without using the leading-log and the mass 
insertion approximations. It proves useful to consider also the "double" ratios, 



R(21/31) 



B(/i^e + 7 ) 
B(r -> e + 7) 



B(r -» e^ r P e ) , R(21/32) = 



B(^^ e + 7 ) 
B(r + 



B(r 



ev T v e 



(17) 



5 Equation 1121 represents the so-called "orthogonal" parametrisation of Y„. In certain cases it is more convenient 
to use the "bi- unitary" parametrisation [23 Y„ = UjjY^ lag Ul, where Ul,r are unitary matrices and YjJ lag is a real 
diagonal matrix. The orthogonal parametrisation is better adapted for our analysis and we will employ it in what 
follows. 
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which are essentially independent of the SUSY parameters. 

To get an estimate for the typical predictions of the schemes with heavy Majorana neutrinos 
with hierarchical spectrum we will consider further, we introduce a "benchmark SUSY scenario" 
defined by the values of the SUSY parameters 



m 



m 



1/2 



250 GeV, A) = a m = -100 GeV , 



(18) 



and tan/3 ~ (5 — 10). In this scenario the lightest supersymmetric particle is a neutralino with 
a mass of ~ 100 GeV. The next to the lightest SUSY particles are the chargino and a second 
neutralino with masses ~ 200 GeV. The squarks have masses in the range of ~ (400 — 600) GeV. 
Supersymmetric particles possessing the indicated masses can be observed in the experiments under 
preparation at the LHC. 

The "benchmark" values of mo, m^j and Aq in eq. (|T8|) correspond to 



B(/i -> + ^ 9.1 x 10 



-10 



y!iy, 



tan 2 (3 , 



(19) 



where (L)^ = ^^(L)^, Lf. = ln(Mx/M^). Since tan 2 j3 will typically enhance B(/j — > lj + 7) 

by at least 1 order of magnitude, the quantity (YtLY v ^j has to be relatively small to be in 

agreement with the existing experimental upper limit on B(^, — > e+7). For given values of the heavy 
Majorana neutrino masses, this will lead to certain constraints on the parameters in R. Regarding 
the masses of the heavy Majorana neutrinos, we shall assume that Mi <C M 2 <C M3, with M3 
having a value M 3 <; (10 13 — 10 14 ) GeV, M3 <C Mx- Constraints from thermal leptogenesis require 
that Mi > 10 9 GeV |HB1 E3J - This would indicate a hierarchy, e.g., of the form Mi ~ (10 9 - 10 
GeV, M 2 ~ (10 12 - 10 13 ) GeV and M 3 ~ (10 14 - 10 15 ) GeV. 



llN 



2.4 Leptogenesis 

In the case of Mi <C M2 <C M3, the baryon asymmetry of interest is given by 

Y B ~ -10- 2 Kei , 



(20) 



where e\ is the CP-violating asymmetry in the decay of the lightest RH Majorana neutrino N\ 
having the mass Mi, and k is an efficiency factor calculated by solving the Boltzmann equations 
(see, e.g., |33j ) . A simple approximate expression for the efficiency factor k in the case of thermal 
leptogenesis we will assume in what follows, was found in |34| : 



1 3.3 x 10" 3 eV 



+ 



mi 



mi 



0.55 x 10" 3 eV 



where the neutrino mass parameter mi is given by 



1.16 



(21) 



mi 



3l (v Vt 

Mi K v 



The CP-violating decay asymmetry ei has the form 



ei 



3 

8^ 



1 



■ V ) 
/ 11 



Im 



Y YtV \ ^ 
YvXv )2i] M 2 



(22) 



(23) 
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Extensive numerical studies have shown joMl IH4| that in MSSM and for hierarchical spectrum of 
masses of the heavy Majorana neutrinos under discussion, successful thermal leptogenesis is possible 
only for 

mi < 0.12 eV . (24) 
For typical values of n ~ one gets for Yb a value compatible with the observations [2] , 

Y B = (6.15 ± 0.25) x 10~ 10 , (25) 

if ei ~ -(10~ 5 - 10~ 7 ). 

As it follows from eqs. ()13|) and ()15|) . the branching ratios of li — > lj + 7 decays in the case 
of interest depend on the orthogonal matrix R. Successful leptogenesis can take place only if R is 
complex, so we will consider (R)* 7^ R. In what follows we will use a parametrisation of R with 
complex angles (see, e.g., |23l l24"|): 



R — R12 R13 R-23 i or R — R12 R23 R- 12 , (26) 

where Rjj (R'i2) describes now the rotation with a complex angle uj-ij = p^ + ioij (u>' 12 = p'\i + 
icr' 12 ), Pij and aij (p' 12 and <t^ 2 ) being real parameters. These parametrisations prove particularly 
convenient for investigating the case of hierarchical spectrum of masses of the heavy RH neutrinos. 

3 The See-Saw Mechanism, Neutrino Yukawa Couplings, LFV Decays 1^ — > Ij+J 
and Leptogenesis 

There has been a considerable theoretical effort in recent years to understand possible connections 
between the neutrino mass and mixing data, LFV charged lepton decays and leptogenesis. Here 
we shall focus on the combined constraints which the existing stringent upper limit on the p — > 
e + 7 decay rate and the requirement of successful leptogenesis impose on MSSMRN. We will be 
interested, in particular, in the possible implications of these constraints for the form of the matrix 
R, the heavy Majorana neutrino masses, the predicted rates of the decays p — > e + 7, r^e + 7 
and r — > p + 7, and the basic SUSY parameters. 

3.1 Normal Hierarchical Light Neutrino Mass Spectrum 

We set first mi, Mi and M 2 to zero. In this approximation we find using R = R12R13R23: 

L3 M3 cos U12 cos 0^3 



(yJly„) 



21 y/2l 



.2 



'm2 I e l ( a ^ M )/ 2 ^/m3 cos 1^23 — \/m 2 C12 sinu>2 3 I S12 S111W23 



-7TT-3 COSW23 COSW23 S13 e l<5 



(27) 



In deriving eq. I|27[) we have set for simplicity #23 = vr/4 and neglected the terms oc m 2 si3 an d 
oc m2 3S?Q in the square brackets. The corresponding expressions for (yJlY v ) are very 

' V 7 31,32 

similar to that for ^YjLY^ . In particular, both are proportional to L3-M3 cosu;i3 cos 0^3. 
For the plausible values of M 3 ^ (10 14 - 10 15 ) GeV and M x = 2 X 10 16 GeV, one finds that 
M3 L3 \ / Am\/(\^2Vu) ~ 0.66. Barring accidental cancellations between the terms in the square 
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least by a factor of ~ 10 3 than the existing upper bound B(/x 
M x = 2 x 10 16 GeV one finds M 3 L 



brackets in eq. (|2*7j) . we get from eq. (|T9*|) that the predicted value of B(/j, — > e + 7) will be larger at 

e + 7). For M 3 ^ 10 13 GeV and 

3 u z y Am|/( v / 2^) = 0.09, which for tan 2 (5 = 25 (100) and the 

chosen "benchmark" values of the SUSY parameters leads to B(/x — > e + 7) = 1.8 x 10~ 10 (7.4 x 
10~ 10 ). The values obtained are still larger than the current limit. This might suggest that the 
SUSY parameter m^j has a bigger value than the "benchmark scenario" value we have assumed, or 
that mo ~ Ti-i/2 <i 500 GeV 6 . We will pursue, however, an alternative hypothesis. We will suppose 



that m 



1/2 



1/2 

few x 100 GeV. If indeed M 3 > 5 x 10 13 GeV, M 3 < M x , where M x > M GU T, the 
existing stringent experimental upper limit on B(/u — > e + 7) might suggest that W13 = tt/2 and we 
will explore this interesting possibility in what follows. For ux3 = tt/2, the R matrix has the form 



sin uj cos u> 
R ~ I cos uj — sin u 
-10 



(28) 



where uj = ojx2 — <^23- Thus, only the combination (oj\2 — 0J23) of the complex angles uix2 and U23 
appears in the expression for R. It is not difficult to convince oneself that if mi is negligible, 
and R has the form given in eq. (|28jl. we have (Y u )^j = (j = 1,2,3). This means that the 
heaviest RH Majorana neutrino N3 decouples and Y v coincides in form with the matrix of neutrino 
Yukawa couplings in the so-called "3x2" see-saw model |48| . We will keep, however, the elements 
(Y l/ ) lj ( 2 j) 7^ (j = 1,2,3) in our further analysis. 

With mi = and R having the form (|28|). the terms ~ M2L2 give the dominant contribution 



111 



yIly l 



1.1 



(i^j). We get: 



YlLY u 



21 



L 2 M 2 



V^Ci 3 S12 C W - e -^M-25)/2 ^ Si3 ^ 



x ( e^l 2 Jm^c 13 S23 s* - e ia/2 Jm^cn c 23 c* 



(29) 



where c w = cosw, s w = sin uj and we have neglected terms oc S13 which can give a correction 
not exceeding approximately 13%. Setting, for instance M2 = 10 12 GeV, we find |(YjLY I/ )2i| ~ 

M 2 L 2 ^ Am\/vl ~ 10~ 2 , and correspondingly B(/z — >■ e + 7) = 2.3 x 1CT 12 for tan 2 /3 = 25. This is 
the range that will be explored by the experiment MEG ^Hj currently under preparation. Similarly, 

we obtain for ( y\,LY v ) : 

V / 31,32 



YILY U 



31 



L 2 M 2 



-ia/2 



'm 2 3X2 c u 



/ ^3C 23 + e 



ia/2 



-i(0M-2S)/2 



>m 2 CX2 S23 c u 



/ m 3 sx3 s L 



(30) 



and 



YlLY, 



■32 



L 2 M 2 



e ia/2 V^2Ci 2 c 23 c w + e iPu/2 Jm^s 2 3S u 



x l e i/W2 c 23 s* + e iQ/2 ^c 12 s 23 < 



(31) 



6 Information about the SUSY parameters mo and mi/j of interest is expected to be obtained in the experiments 
under preparation at the LHC. 
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As it follows from eqs. (|29j ) -l|31j ) . for u> / 0, the k — > lj +7 decay branching ratios of interest depend 
on the Majorana phase difference (a — Pm)- The effective Majorana mass \{m)\ in (/3/3)oj,-decay 
(see, e.g., ^UEI]) depends on the same Majorana phase difference (see, e.g., [3H1IBH]) : 



|{m> I 



Am\\ sin 2 9i2e 



i(a-j3 M ) 



+ 



im^ sin 2 



(32) 



If s^ls^l is not negligibly small, B(p — > e + 7) and B(r — > e + 7) will depend also on the Dirac 
phase 5. 

For the "double" ratio R(21/31) we find from eqs. (l29l) and (pp : 



R(21/31) 



/m 3 s 2 3 s u 



/ rn2C 12 c 23 c UJ e l 2 



/mi c 23 + Jm2~ci2 s 2 3 c w e J 



(33) 



Given #12, #23 arid m2/m 3 = ^/Amfi/yAm^, R(21/31) depends only on (a — /3m) and w. If the 
terms oc ^/to3 (oc y / mi) in eq. (fBl^ dominate, we have R(21/31) = 1. 

The expression for the double ratio R(21/32) can be obtained from eqs. (|29j) and (|3*T1) . For 
\frnisi3,\su\ > we get R(21/32) = s^/c^ Si 0.1, while if ^/mi\sj\ <C v 7 "^"! I s i2 , on e 

finds R(21/32) = (tan 2 6>i 2 )/s| 3 0.9. 



3.1.1 Leptogenesis Constraints 

We shall analyse next the constraints on the parameter oj in the matrix R, eq. I)28|) . which follow 
from the requirement of successful thermal leptogenesis. With to\ 3 = tt/2 we find that in the case 
of NH light neutrino mass spectrum we are considering. 



3 / m 3 M 1 



Im 



2 Am|i 2 



Am 



31 



|cj 2 + 



3 

8^ 



m 3 M x 



sin 2p sinh2<T 



(1 



^)cos2p+(l + ^)cosh2a 



(34) 



(35) 



where /) and a are determined by u = p + id and we have used m\ = Am^, m\ = Am^ and the 
relation Imc 2 = — Ims 2 . Thus, in the case under discussion, the mass M2 governs the magnitude of 
the li — ► lj + 7 decay branching ratios, whereas Mi determines the value of the leptogenesis decay 
asymmetry. The conditions | | 2 > 0, | | 2 > imply that p and a (in eq. l|3*5)) ) should satisfy 
cosh2cr > I cos2/)|, which is always valid since cosh2cr > 1. As can be easily shown, we have 



Im 



c?„ + 



Am 21 ,2 



Am 



31 



< 



Imc 



< 1 



cJ + 



, ,2 
AmJ, ^ 



(36) 



leading to the well-known [IS] upper limit 

3 /m 3 M : 



lea I ^ 



87T \ U 2 



* ~ 1.97 x 10~ 7 



m 3 



0.05 eV/ \ 10 9 GeV 



Mi 



174 GeVY 



(37) 
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Figure 1: The leptogenesis constraints on the parameters p and a for Mi = 10 10 GeV; 10 11 GeV; 
10 12 GeV (blue; green; red areas) in the case of NH light neutrino mass spectrum. The two 
panels correspond to two different intervals of values of the baryon asymmetry of the Universe, Yb, 
considered: (a) 5.0 x 10" 10 <Y B < 7.0 x 10~ 10 , and (b) 3.0 x 10~ 10 <Y B < 9.0 x 10~ 10 . The solid 
lines show the limit associated with the upper bound rh\ < 0.12 eV: outside the region between 
the solid lines the wash-out effects are too strong and leptogenesis cannot produce the observed 
baryon asymmetry. 



The requirement of a nonzero asymmetry, ei 7^ 0, implies, as it follows from eq. (|35[). that both the 
real and imaginary parts of u have to be nonzero, p 7^ kir/2, k = 0, 1, 2, ff^O, i.e., that R has 
to be complex. Moreover, we should have sin 2psinh 2a < since the decay asymmetry e\ has to be 
negative in order to generate a baryon asymmetry of the correct sign. The maximal asymmetry |ei| 
is obtained for |Imc 2 | = |c w | 2 , which is satisfied for cos2p cosh2<7 = — 1, cos2p 7^ — 1, cosh2<r 7^ 1. 
The neutrino mass parameter fh\, eq. (|22|). can also be easily found: 

1 1 

mi ~ m 3 |c^| 2 + m2\s LU \ 2 = ^(""^ + m 2) cosh2a + ^("^ ~ ^2) cos 2p > m 2 . (38) 

The minimal value of fh\ = m,2 = yAmfi = 9 x 10 -3 eV, corresponds to cosh2cr = 1 and 
cos2p = —1, for which |ei| = 0. For 9 x 10~ 3 eV < rh\ ^ 0.12 eV, where we have taken into 
account eq. (|24|). the efficiency factor lies in the interval 1.9 x 10~ 3 5s k < 3.9 x 10~ 2 . For this 
range of values of k successful leptogenesis is possible for Mi ^ 10 10 GeV. We will consider values 
of Mi in the interval Mi = (10 10 — 10 12 ) GeV, which is compatible with the assumption we made 
about the hierarchical mass spectrum of the heavy Majorana neutrinos. Thus, for given Mi, the 
requirement of successful leptogenesis implies a constraint on the two parameters p and a of the 
theory. In Fig. ^ we show the leptogenesis constraint on p and a for Mi = 10 10 GeV; 10 11 GeV; 
10 12 GeV. As we see from Fig. ^ the requirement of successful leptogenesis severely limits the 
allowed ranges of values of p and a. Moreover, the values of the two parameters are strongly 
correlated. We note, in particular, that as |<r| increases, the wash-out effects become stronger and 
for \a\ ^ 1, the observed baryon asymmetry cannot be reproduced. The maximal asymmetry |ei| 
for a given Mi is obtained for values of cos 2p and cosh2<r close, but not equal, to (—1) and (+1), 
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Figure 2: The correlation between the predicted values of the LFV decay branching ratios B(/j — > 
e + 7) (a), B(r — > e + 7) (b), B(r — > \i + 7) (c) and of the baryon asymmetry Yb in the thermal 
leptogenesis scenario, for M\ = 6x 10 10 GeV, M2 = 10 12 GeV and NH light neutrino mass spectrum. 
The figure was obtained for the "benchmark" values of the soft SUSY breaking parameters ttiq = 
m i/2 = 250 GeV, CLQrriQ = —100 GeV and the minimal value of tan/3 = 5. The region between 
the two vertical dashed lines corresponds to the observed baryon asymmetry: 5.0 x 10~ 10 < Yb < 
7.0 x 10~ 10 . Results for two values of the Majorana phase (a — /3m) equal to (red+green areas) 
and 7r (green areas) are shown. 



respectively. Given the interval of allowed values of mi, we can write rh\ = f 771,3 with / = [0.2,2]. 
The condition of maximal |ei| implies cos 2/7 = (/ — \J\ + / 2 ) 7773/ (777.3 — m 2)- Choosing / = 1, i.e., 
777-1 = m 3 — 0.05 eV, for instance, we get cos2p = —0.5 and correspondingly cosh2er = 2. In this 
example k = 5.4 x 10 -3 and for tan 2 (3 <; 3 we get the requisite value of the the baryon asymmetry 
for Mi 6 x 10 10 GeV. 

In Fig. |21 we show the relation between the predicted values of Yb in the thermal leptogenesis 
scenario and of B(Zj — > The figure was obtained for the "benchmark" values of the soft SUSY 

breaking parameters ttiq = 777 1 / 2 = 250 GeV, ao m o — —100 GeV and the minimal value of tan/3 = 5. 
For the relevant heavy Majorana neutrino masses we used Mi = 6 x 10 10 GeV and M2 = 10 12 GeV. 
Results for two values of the Majorana phase difference {a — (3m) = 0; it, are shown. For the values 
of p and a ensuring successful leptogenesis we find that typically 10 -14 ^ B(/i — > e + 7) ^ 5x 10~ 13 , 
10- 13 < B(t -> n + 7) < 5 x 10~ 12 and lO" 15 < B(r e + 7) < 5 x 10~ 13 . However, we have 
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B(li — > /j + 7) oc tan 2 /3 and, e.g., for tan/3 = 20 we get typically 1.6 x 10 -13 Ss B(/z — ► e + 7) 
8 x 10~ 12 , which is entirely in the range of sensitivity of the MEG experiment. As Fig. ^indicates, 
the dependence of B(Zj — » Zj + 7) on the Majorana phase (a — (3m) is relatively weak. 

3.2 Inverted Hierarchical Light Neutrino Mass Spectrum 

We will perform next a similar analysis assuming that the light neutrino mass spectrum is of the 
inverted hierarchical type. We set m-2 = mi, m^/mx^ — and neglect first M1/M3 and M2/M3. 
Setting for simplicity S13 = and #23 = 7r/4, we find that B(/i — ► e + 7) will depend on 

/ t \ -L 3 M 3 m 2 ( . _ ia/2 . \ 

[YILY U J^ ~ — \ci2 smaji3 + si 2 e 1 cos w 13 sm u; 23 J 

x [e ia l 2 c\2 cosw* 3 sin 1^23 — «12 si nCt; i3^ • (39) 
This serves to underline that - as in the case of a NH light neutrino spectrum discussed in Section l3~Tl 



we have typically 



- M 3 ^Am 2 A /vl. This leads for M 3 ~ (10 14 - 10 15 ) GeV 

and mo, Aq and mi^ in the fewxlOO GeV range to a n — > e + 7 decay branching ratio which 
exceeds the existing limit by approximately 3 orders of magnitude. Looking again for simplifications 
with interesting phenomenological consequences, we can reduce the magnitude of (yIlyJ^ by 

setting W13 = u>23 = and, correspondingly, R13 = R23 = 1 in eq. (|2*o]) 7 . The corresponding form 
of R is 

(COSW12 sinu;i2 \ 
— sinu;i2 COSW12 . (40) 
01/ 

With R given by eq. (|40[) and negligible m 3 /mi2, the heaviest (RH) Majorana neutrino .ZV3 decou- 
ples and we have again (Y u ) 3 j = (j = 1,2,3). Neglecting further the splitting between mi and 
m.2 we find: 



/ vtrv \ L 2 M2^J\Aml 1 \ , , 
[YlLYyj^ ~ ^ \ci2 smu\2 - e 1 s 12 cosw^ 



e ta/2 C12 C23 cos u*2 + (S12 C23 + e lS C12 S23 s 13 ) sina;* 2 



(41) 



where we have used m.^2 — ^/(Am^-J and have neglected terms oc S13 which give a correction 



not bigger than approximately 13%. Being of the order M2\ \ Am^ \/v 2 , the expression (|4*T|) for 



YjLY^ ) will lead for M 2 ~ 10 12 GeV and values of the soft SUSY breaking parameters in the 
/ 21 

fewxlOO GeV range to B(u — > e + 7) close to the existing limits. For ( YiLY u ) we similarly 

V / 31,32 



7 If one uses a somewhat different parametrisation of R, namely, R = Ri2R23R'i2, the same result is achieved by 
setting just tJ23 = 0. 
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get 



31 



L 2 M 2 J Amf! 



-ia/2 



S12 cosa;i2 - c i2 sina)i 2 



e 7 C12 S23 COS(J 12 + (Si2 S23 - e C12 c 23 S13) smw 



12 



(42) 



32 



L2M2A/ Amfj 



C23 S12 + e * 5 C12 S23 S13 ) sin a; 12 + e 40/2 c i2 c 23 cos uj 12 



e 7 C12 S23 cos a; 12 + si2«23-e C12 c 23 si 3 sinw 12 



(43) 



We see that, as in the case of NH light neutrino mass spectrum, ^YjLY^j and ^YjLY^j are 

rather similar in structure, whereas ( yJlY,, ) differs somewhat. The Majorana phase Bm does 

V / 32 

not appear in the expressions (|41j) - (|43|) because we have set 777,3/771,^2 = 0. As the phase factor 
including the Dirac phase 5 appears always multiplied by the small parameter S13, for S13 < 0.1 
the branching ratios depend essentially only on the Majorana phase a, which enters also into the 
expression for the effective Majorana mass |(t7i)| in (/3/3)oj,-decay [5U1 1351 13*§] : 



Am 2 3 |cos 2 0i2 + e ia sin 2 1 



'12 



(44) 



We will give next the ratios of B(/x — > e + 7) and B(r — > e + 7) (B(r — > /_i + 7)) in the case of 
negligible contribution of the terms oc S13 8 : 



R(21/31) ~ cot 2 (9 23 



R(21/32) ~ s 



-2 

23 



Ja/2 



C12 smtJi2 — S12 cosa;i2 



?f> /2 

C12 sinwi2 - e 7 S12 cosa;i2 



3 ia/2 



S12 sino;i2 + ci 2 COSW12 



?7> /2 

si2sm(Ji2 + e 7 C12 COSW12 



(45) 



(46) 



Hence, as in the case of NH light neutrino mass spectrum, R(21/31) is rather close to one, whereas 
R(21/32) can have a wide range of values. Most interestingly, R(21/32) can have a value close to 
two or even be as large as ~ 10. 



3.2.1 Leptogenesis Constraints 

We can again work out possible constraints from the requirement of successful leptogenesis. Using 
expression (j40jl for the matrix R and eq. (|23|) we find that the CP-violating decay asymmetry e% 
of interest has the form 



3 f m 2 Mi\ Ami 
61 - " — ( — — 



8tt 



Im [sin 2 W12] 



„ / lAmfJ 



n Am^ \ , . |2| , ,2 
1 + I . 2 I I smw i2| +|cosa;i2| 
2 Am 31 / 



3 / / m 2 Mi^ Am^i 
16vr 



Am§J 



sin2p tanh2cr , 



(47) 



3 For S13 < 0.10 the correction due to the terms in question can be shown to be smaller than approximately 15%. 
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Figure 3: The same as in Fig. ^ but for IH light neutrino mass spectrum and (a) M\ = 7 x 10 12 
GeV; 1.5 x 10 13 GeV; 3 x 10 13 GeV (blue; green; red areas), (b) M x = 4 x 10 12 GeV; 7.0 x 10 12 
GeV; 1.5 x 10 13 GeV (blue; green; red areas). 



where oj\2 = p + icr, mi = y^Am 2 ^ and we have neglected corrections ~ Am^/ lAmlJ. We see 
that in order to have ei 7^ 0, both p and a should be different from zero: p 7^ kn/2, k = 0, 1, 2..., 
a 7^ 0. Since e\ < 0, we should have sin 2/otanh 2a > 0. 

It follows from eq. (|47j) that in the case of IH light neutrino mass spectrum under discussion, the 
CP-asymmetry t\ is suppressed by the factor Ara^/ | Am^ | . The expression for the CP-asymmetry 
we have found for the NH spectrum, eq. 1)35 jl . does not contain the indicated suppression factor. It 
is not difficult to show that one always has 

|lm[sin 2 u, 12 ]| ^1 (4g) 



Am 2 \ , |2 , , ,2 ~ 

1 + ~Ta — 21 p mi ^i2 + cosu;i2 



For the asymmetry ei we get the upper limit 

, ei | < J_ (H^l) „ 3 2 x 10- 9 f^M ( Ml ) ( 1U GeV V (49) 

|ei| ~16vrV vl ) \Aml\~ V0.05 eW VlO 9 GeV/ V v u ) ' { ^ 

where we have used Am^/ lAm^j = 3.2 x 10~ 2 . The maximal value of ei is reached for p = tt/4 
and a <^ 0.5. 

For the neutrino mass parameter mi, eq. (|22j) . we find: 

rh\ ~ mi )2 (| cos W12I 2 + | sinu;i2| 2 ) = cosh 2a > mi )2 • (50) 



The minimal value of m\ = = y [Am 2 ^! = 5 x 10 2 eV, corresponds to cosh2<7 = 1 

for which |ei| = 0. For 5 x 10~ 2 eV < rh\ ^0.1 eV, the efficiency factor lies in the interval 
k < 5.4 x 10~ 3 . It is not difficult to convince oneself that for a given Mi, the 
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Figure 4: The correlation between the predicted baryon asymmetry of the Universe Yg and the 
predicted branching ratio of \i — > e + 7 decay B(^ — > e + 7) in the case of IH light neutrino mass 
spectrum and for Mi = 7.0 x 10 12 GeV and M 2 = 4.0 x 10 13 GeV. The fi gure was obtained for the 
"benchmark" values of the soft SUSY breaking parameters mo = mi/2 = 250 GeV, aomo = —100 
GeV and tan/3 = 5. The horizontal line indicates the experimental upper limit on B(^ — ► e + 7), 
while the region between the two vertical dashed lines is favoured by the observed value of the 
baryon asymmetry of the Universe, 5.0 x 10 -10 < Yb < 7.0 x 10 -10 . 



maximal value of |ex| is reached for a = ±0.5, for which rh\ = 7.8 x 10~ 2 eV and, correspondingly, 
k = 3.2 x 10~ 3 . Thus, successful leptogenesis can take place for Mi <; 6.7 x 10 12 GeV, where we 
have used eq. (|25[). In Fig. |3]we show the regions of values of p and a, favoured by requirement 
of successful thermal leptogenesis, for three fixed values of Mi = 7 x 10 12 GeV; 1.5 x 10 13 GeV; 
3 x 10 13 GeV. We find that |<r| ^ 0.75, |<r| 7^ 0. Given the minimal value of Mi determined by the 
leptogenesis constraint, we will consider further in this subsection the following hypothetical heavy 
Majorana neutrino mass spectrum: Mi = 7. x 10 12 GeV, M 2 = 4.0 x 10 13 GeV, M 3 = 2.0 x 10 14 



GeV. For M 2 = 4.0 x 10 13 GeV, we find M 2 J \ Ara^ | L 2 /v^ = 0.4. At the same time the upper limit 

on B(;i ^ e + 7) implies that for the "benchmark" values of the soft SUSY breaking parameters 

2 



we have specified earlier and tan/3 = 5, we should have 



yIly l 



< 4.8 x 10" 4 . It follows 



rt ^ 2 



21 

from the explicit expression for I (YjLY^J |~, eq. (jUJ, that this upper limit is impossible to 

satisfy for the values of the parameters p and a satisfying the leptogenesis constraint (Fig. [2J). 
This is clearly seen in Fig. which shows that the requirement of successful leptogenesis and the 
existing experimental upper limit on B(/x — > e + 7) are incompatible in the case of the "benchmark" 
values of the SUSY parameters, tjiq = m 1 / 2 = 250 GeV, aomo = —100 GeV, and of tan/3 = 5. 
The result we have obtained indicates that in the case of IH light neutrino mass spectrum, the 
SUSY parameters tuq and/or m 1( / 2 should have values considerably larger than the "benchmark" 
values we consider. More specifically, we can have mo ~ (250 — 300) GeV, but <C m 2 ^ 2 . This 
possibility is illustrated in Fig. [51 where the predicted values of B(Zj — > lj + 7) for mo = 300 
GeV, m!/ 2 = 1400 GeV, aomo = 0, and tan/3 = 5, are shown as functions of the predicted value 
of the baryon asymmetry. The figure corresponds to Mi = 7 x 10 12 GeV and M 2 = 4.0 x 10 13 
GeV. Now the requirement for successful leptogenesis is compatible with the existing constraint 
on B(/i — > e + 7): for the values of p and a ensuring successful leptogenesis we find that typically 
3 x 10~ 14 ^ B(p -> e + 7) ^ 5 x 10~ 13 . Significantly larger values of B(p — > e + 7) are possible 
if tan (3 ^ 10. As Fig. also shows, the predicted branching ratios B(Zj — ► L + 7) exhibit weak 
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Figure 5: The correlation between the predicted Yg and the predicted B(/i — ► e+7) (a), B(r — * e+7) 
(b), B(r — > /1+7) (c), for IH light neutrino mass spectrum and Mi = 7.0x 10 13 GeV, M 2 = 4. Ox 10 13 
GeV. Th results shown are obtained for two values of the Majorana phase a = 0; ir and the following 
set of values of the SUSY parameters: niQ = 300 GeV, m 1 / 2 = 1400 GeV, ao = and tan/3 = 5. 
The region between the two vertical dashed lines corresponds to 5.0 x 10~ 10 < Yb < 7.0 x 10~ 10 
and is favoured by the measured value of the baryon asymmetry. The green (red+green) areas 
correspond to the Majorana phase a = n (0). 



dependence on the Majorana phase a. 

If the light neutrino mass spectrum is of the IH type, the results we have obtained in this 
subsection can have important implications for the predicted spectrum of SUSY particles in the 
few 100 GeV - 1 TeV region, to be probed by the experiments at the LHC. For mo = 300 GeV, 
m i/2 = 1400 GeV, aomo = 0, and tan/3 = 5, the lightest SUSY particle is still a neutralino and its 
mass is approximately 600 GeV. The mass of next to the lightest SUSY particle, which is a stau, 
is very close to the mass of the lightest neutralino. At the same time the squarks are predicted to 
be relatively heavy, having masses ~ (2 — 3) TeV. 

3.3 Quasi-Degenerate Light Neutrinos 

In this case one has mi = m 2 = 711,3 = m ) with m ^ 0.1 eV. It is easy to see that ^YjLY^ will 
be proportional to M^m/v^, and therefore a too large branching ratio for [i — ► e + 7 decay will be 
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predicted. Indeed, setting Mi = and M2 = and using the complex Euler angle parametrisation 
R = R' 12 R.23Ri2, we find: 



yIly, 



21 



L 3 mM 3 
V2vl 



(si2 + e l<5 si3Ci2) sinu;^ sino;* 2 — e l 2 C12 coso^ sinw^ + e l a" C13 cos a; 



23 



Ci3(ci2 sinwi2 - e l 2 si 2 coswi 2 ) sina;23 + e" 



S13 COS W 2 3 



(51) 



where for simplicity we have set #23 = tt/4 and have neglected the sub-dominant terms oc S13. 



There are two possibilities for suppression of 

rt 



yJly, 
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: (£) If sina;23 = 0, the contribution 



due to M 3 in I Y V LY v I remains, but is proportional to s\ 3 . Xhe necessary suppression can take 
place if S13 is sufficiently small, (ii) The parameters W12 and CJ23 can have values such that the 
different terms oc M 3 in (yIlyJ^ cancel (completely or partially) each other. The latter seems 
to require fine tuning between the values of several very different parameters. 

In what follows we shall consider the case (i) and we set W23 = 0. In this case R has the form 
given in eq. (|4U|). The quantities of interest ( YjLY^ ) (i 7^ j), including the contributions oc M2, 

\ ) ij 

are given by: 



YtLY v 
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L 3 mM 3 iS 

-e S23C13S13 



i-'r. 



-C13 



^C23S12 + e % S23C12S13)S^ + 2 C23C12C* (— Cl2S w + e 



S12Cuj) , 

(52) 



yJxy, 
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L 3 mM 3 iS 

5 — e C23C13S13 



L2mM2 



-C13 



-S23S12 + e C 23 Ci2Sl3)^ - e 2 S 2 3Cl2C a 



(53) 



Y.tLY, 
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L 3 mM 3 2 

2 c 13 c 23«23 



u 

L2mM2 



-C13 



-S23S12 + e C 23 Ci2Sl3)^ - e 2 S 2 3Cl2C a 



(C235l2 + e * <S S23Cl2Sl3)So; + e i 2C 2 3Ci 2 Q i 



(54) 



where u = uji2+uj' 12 - It is interesting to note that the quantity \Y\,LY v j , and correspondingly 

B(r — ► \i + 7), is not suppressed by the factor sf 3 . The effective Majorana mass in (/3/3)o^-decay 
depends in the case of QD spectrum on the CP- violation Majorana phase a |38[ I39| . present in the 
expressions for (YiLY, 



1.1 



(i^j): \(m)\ ^m\cos 2 e 12 + e ia sin 2 9 12 \ 



3.3.1 Leptogenesis Constraints 

For the CP-violating decay asymmetry ei we find 
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87T \ u 2 / m 
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(55) 
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Figure 6: The same as in Fig. ^ but for QD light neutrino mass spectrum and (a) M\ = 3.0 x 10 13 
GeV; 5.0 x 10 13 GeV (green; red areas), (b) M x = 1.7 x 10 13 GeV; 3.0 x 10 13 GeV (green; red 
areas) . 



It is not difficult to show that 

Im(4) 



Im( S 2) 



M 2 + (i + ^)k,| 2 + k 



T7 = — sin 2/9 tanh 2cr < — 
I 2 2 ~ 2 



(56) 



Thus, the maximal asymmetry |ei| is given by 

'0.1 eV 



ei <1.6 x 10 



-9 



Mi 



174 GeV 



to j \ 10 9 GeV . 
One can easily find also the mass parameter rh\. 

m, cosh 2c > to . 



TOl = TO 



I I 2 _i_ I |2 

t Sin 



(57) 



(58) 



Since successful leptogenesis is possible for |34j fh\ < 0.12 eV, while for QD light neutrino mass 
spectrum to ^ 0.1 eV, we get from eq. ()58(l that to = 0.1 eV. Therefore in all further analysis and 
numerical calculations in this subsection we set m = 0.1 eV. 

As it follows from the preceding discussion, we have rh\ = (0.10 — 0.12) eV. Correspondingly, the 
wash-out effect in the case under consideration is relatively strong. Taking into account the precise 
upper limit on rh\ given in ref. [34] : fh\ < 0.12 eV, we get for the corresponding efficiency factor 
1.9 x 10 -3 <k<2.4x 10 -3 . The condition rh\ < 0.12 eV implies a ^ 0.3, for which tanh2cr < 0.5. 
Thus, using eqs. (|2T?)) and (|57j) we obtain the minimal value of M\ ensuring successful leptogenesis: 
M\ 3.0 x 10 13 GeV. Since rh\ ft 0.1 eV, we can expect that a lies in the interval a = (0.2 — 0.3). 
This is confirmed by a more detailed numerical analysis. The values of the parameters p and a 
allowed by the leptogenesis constraint are shown in Fig. ©for Mi = 3.0 x 10 13 GeV; 5.0 x 10 13 GeV. 

Given the lower bound M\ <; 3.0 x 10 13 GeV, a possible mildly hierarchical heavy Majorana 
neutrino mass spectrum would correspond to, e.g., M\ = 3.0 x 10 13 GeV, M2 = 1.2 x 10 14 GeV, 
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Figure 7: The correlation between the predicted Yq and the predicted B(/j, — > e + 7) in the 
case of QD light neutrino mass spectrum and for Mi = 3.0 x 10 13 GeV, M 2 = 1.2 x 10 14 GeV 
and M3 = 4.8 x 10 14 GeV. The "benchmark" values of the soft SUSY breaking parameters mo = 
m.1/2 = 250 GeV, clqttiq = —100 GeV and tan/3 = 5, have been used. The horizontal line indicates 
the experimental upper limit on B(/j — > e + 7), while the region between the two vertical dashed 
lines corresponds to 5.0 x 10~ 10 < Y# < 7.0 x 10~ 10 and is favoured by the observed value of Yg. 



and M 3 = 4.8 x 10 14 GeV. For this spectrum, L 3 mM 3 /vl = 6.0 and L 2 mM 2 /vl = 2.0. Using 
the lowest possible value for tan 2 j3 = 10 we find that even if S13 = and the term oc M3 does 



not contribute to 



yUy v 



the contribution of the terms oc M 2 is so large in the case of the 

'21 

"benchmark" values of the soft SUSY breaking parameters, mo = 7771/2 = 250 GeV, ao777o = —100 
GeV, that the predicted B(/i — ► e + 7) exceeds the existing upper limit on B(/x — ► e + 7) by more 
than a factor of ~ 10 3 . The indicated incompatibility between the leptogenesis and B(/i — ► e + 7) 
constraints is illustrated in Fig. [7J 

Similarly to the case of IH light neutrino mass spectrum, the requirement of successful thermal 
leptogenesis and the upper limit on B(/i — > e + 7) can be simultaneously satisfied only if the scale of 
masses of supersymmetric particles is significantly larger than that predicted for the "benchmark" 
values of the soft SUSY breaking parameters we have adopted. In Fig. |H] we show the correlation 
between the predicted values of B(Zj — ► lj + 7) for sin 6*13 = 0.05 and tan/3 = 5, mo = 300 
GeV, mi/2 = 1400 GeV and ao = 0, and the predicted value of the baryon asymmetry. The results 
presented in this figure have been obtained for the spectrum of the heavy Major ana neutrino masses 
specified above. We note, in particular, that the predicted interval of values of B(/i — ► e + 7) which 
is compatible with the observed baryon asymmetry is in the range of sensitivity of the ongoing MEG 
experiment: B(/i — > e + 7) can have a value just below the present experimental upper limit. As in 
the cases of NH and IH light neutrino mass spectra, we find that the dependence of — ► lj + 7) 
on the relevant Major ana phase a is rather weak. 

In Fig. the correlations between the predicted value of Yb and those of the double ratios 
R(21/31) and R(21/32) are displayed. In the approximation we use, the double ratios are in- 
dependent of SUSY parameters and are determined only by the off-diagonal elements of Y\Y V . 
When the constraint of successful leptogenesis is imposed, we get for the allowed range of values of 
R(21/31) for the NH, IH and QD light neutrino mass spectrum respectively 10 ^ R(21/31) ^ 100, 
20 < R(21/31) ^ 50 and R(21/31) ~ 20; 100. Similarly, for the double ratio R(21/32) we get in the 
three cases lO" 2 < R(21/32) < 10, 10 < R(21/32) < 10 3 and 50 < R(21/32) < 10 3 , respectively. 
We find, in particular, that R(21/32) can be much smaller than 1 only for NH light neutrino mass 
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Figure 8: The correlation between the predicted Yb and the predicted B(/i — > e + 7) (a), B(r —* 
e + 7) (b), B(r — > \i + 7) (c), for QD light neutrino mass spectrum and Mi = 3.0 x 10 13 GeV, 
M 2 = 1.2 x 10 14 GeV, M 2 = 4.8 x 10 14 GeV. The SUSY parameters used to obtain the figure 
are too = 300 GeV, m 1 / 2 = 1400 GeV, ao = and tan/3 = 5. Results for two values of the 
Majorana phase a = 0; ir are shown. The region between the two vertical dashed lines corresponds 
to 5.0 x 10~ 10 <Y B < 7.0 x 1Q- 10 . 



spectrum. 

4 Conclusions 

We have considered the LFV decays /i^e + 7, r^e + 7 and r — ► fj, + 7 and leptogenesis 
in the MSSM with see-saw mechanism of neutrino mass generation and soft SUSY breaking with 
universal boundary conditions at a scale Mx > Mr, Mr being the heavy (RH) Majorana neutrino 
mass scale. The heavy Majorana neutrinos were assumed to have hierarchical mass spectrum, 
Mi < M 2 <C M 3 , while the scale M x was taken to be the GUT scale, M x = 2 x 10 16 GeV. We have 
derived the combined constraints, which the existing stringent upper limit on the fi — > e + 7 decay 
rate and the requirement of successful leptogenesis impose on the neutrino Yukawa couplings, heavy 
Majorana neutrino masses and SUSY parameters in the cases of the three types of light neutrino 
mass spectrum - normal and inverted hierarchical (NH and IH), and quasi-degenerate (QD). A 
basic quantity in these analyses is the matrix of neutrino Yukawa couplings, Y u . In the present 
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Figure 9: The correlation between the predicted Yb and the double ratios R(21/31) and R(21/32), 
defined in eq. (|17|) . The heavy neutrino mass spectrum in the NH, IH, and QD cases is taken to 
be same as in Fig. |21 Fig. El and Fig. |H1 respectively. The region between the two vertical dashed 
lines corresponds to 5.0 x 1CT 10 <Y B < 7.0 x 10~ 10 . 



work we have used the orthogonal parametrisation of Y„, in which ~Y V is expressed in terms of the 
light neutrino and heavy RH neutrino masses, the PMNS neutrino mixing matrix ?7pmns> an d an 
orthogonal matrix R. Leptogenesis can take place only if R is complex. 

The constraints from thermal leptogenesis require, in general, that M x > 10 9 GeV. This would 
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indicate a hierarchy, e.g., of the form Mi ~ (10 9 - 10 11 ) GeV, M 2 ~ (10 12 - 10 13 ) GeV and 
M 3 <; 10 13 GeV » M 2 , M 3 < (<)Mjf. In our analysis we have considered a "benchmark SUSY 
scenario" defined by the values of the soft SUSY breaking parameters in the range of few x 100 
GeV: mo = = 250 GeV, aQiriQ = —100 GeV, and tan/3 ~ (5 — 10). In this scenario the 

lightest supersymmetric particle is a neutralino with a mass of ~ 100 GeV. The next to the lightest 
SUSY particles are the chargino and a second neutralino with masses ~ 200 GeV. The squarks 
have masses in the range of ~ (400 — 600) GeV. Using the indicated set of "benchmark" values 
of the soft SUSY breaking parameters and barring accidental cancellations, we find that for the 
typical values of the heaviest Majorana neutrino mass M3 = (5 x 10 13 — 10 15 ) GeV, the limit on 
the n — > e + 7 decay branching ratio B(/u — > e + 7) is impossible to respect independently of the 
type of the light neutrino mass spectrum: the predicted B(/j, — > e + 7) exceeds the existing upper 
limit by few orders of magnitude. For each of the three types of neutrino mass spectrum - NH, 
IH and QD, we find simple forms of the matrix R which lead to a suppression of the dominant 
contributions due to the terms oc M3 in B(/i — ► e + 7). In all three cases the matrix R ensuring the 
requisite suppression admits a parametrisation by one complex angle. In the case of NH spectrum 
R is given by eq. (|28|) . while for IH and QD spectra, R = R12, R12 being the matrix of (complex) 
rotations in the 1-2 plane. In this case the dominant contribution in B(// — > e + 7) comes from 
terms oc M 2 . For QD spectrum the terms oc M3 are suppressed by the factor sin #13 and can be 
comparable to those oc M 2 . 

The requirement of successful leptogenesis leads to a rather stringent constraint on the complex 
mixing angle in R. For IH and QD spectra it also implies a relatively large lower limit on the 
mass of the lightest RH Majorana neutrino: Ml > 7.0 x 10 12 GeV and M x > 3.0 x 10 13 GeV, 
respectively. With such values of Mi and hierarchical heavy Majorana neutrino mass spectrum, 
the upper bound on B(^ -> e + 7) can be satisfied only if the scale of masses of SUSY particles 
is considerably higher than that implied by the "benchmark" values of the soft SUSY breaking 
parameters we have considered. We have analysed a specific case of such SUSY scenario: uiq = 300 
GeV, mi/2 = 1400 GeV and ao = 0. In this scenario the lightest SUSY particle is a neutralino with 
a mass of approximately 600 GeV, the next to the lightest SUSY particle is a stau and its mass 
is very close to the mass of the lightest neutralino, while the squarks are relatively heavy, having 
masses ~ (2 — 3) TeV. The predictions for B(/i — > e + 7) are now largely in the range of sensitivity 
of the ongoing MEG experiment. If more stringent upper limits on B(/x — > e + 7) will be obtained 
in the future, it would be rather difficult to reconcile the IH and QD light neutrino mass spectra 
with the n — > e + 7 and leptogenesis constraints and SUSY particle masses in the TeV range. Our 
results may have important implications for the search of SUSY particles in the fewxlOO GeV - 1 
TeV region, to be performed by the experiments at the LHC. 

Satisfying the combined constraints from the existing upper limit on the /i — > e + 7 decay rate 
and the requirement of successful thermal leptogenesis proves to be a powerful tool to test the 
viability of supersymmetric theories with see-saw mechanism of neutrino mass generation and soft 
flavour-universal SUSY breaking at a scale above the heavy RH Majorana neutrino mass scale. 
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